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FRACTAL ANALYSIS OF HOPF BIFURCATION AT 

INFINITY 

GORAN RADUNOVIC, DARKO ZUBRINIC AND VESNA ZUPANOVIC 


Abstract. Using geometric inversion with respect to the origin we ex¬ 
tend the definition of box dimension to the case of unbounded subsets of 
Euclidean spaces. Alternative but equivalent definition is provided us¬ 
ing stereographic projection on the Riemann sphere. We study its basic 
properties, and apply it to the study of the Hopf-Takens bifurcation at 
infinity. 


1. Introduction 


The main goal of dimension theory for dynamics is to measure the com¬ 
plexity of invariant sets and measures using fractal dimensions. A good 
example of this can be seen in Piacquadio et a/. (2002)] where the Hausdorff 
dimension of a particular case of the Henon attractor is estimated and com¬ 
pared to its box dimension. In many cases fractal dimensions can give us 
a better understanding of the dynamics appearing in various problems in 
physics, engineering, chemistry, medicine, etc. The second and the third au¬ 
thor used th e box dimension to analyse spir al trajectories of some planar vec¬ 


tor fields in Zubrinic &: Zupanovic(2005) . Among other things, they stud¬ 


ied the Hopf bifurcation which is a well known bifurcation of 1-parameter 
families of vector fields in which a limit cycle, that is, an isolated periodic 
orbit, is born from a singular point. The generalisation is the Hopf-Takens 
bifurcation which gives rise to more than one limit cycle born from a singular 
point. It was shown that the box dimension of spiral trajectories near singu¬ 
lar points or nonhyperbolic (multiple) limit cycles becomes nontrivial, that 
is, greater than 1 precisely at a point at which the corresponding dynamical 
system undergoes the bifuraction. Moreover, the box dimension can only 
take values from a discrete set and depends on the multiplicity of the corre¬ 
sponding singular point or limit cycle. This could be utilised for computing 
the multiplicity. In a way, this is related to the 16th Hilbert problem of hnd- 
ing an upper uniform bound for the number of limit cycles in dependence on 
the degree of the polynomial vector held. R esults about spiral trajectorie s 
of some vector helds in can be found in [Zubrinic &: Zupanovic(2005) 


Further studies using the asympthotic behaviour of the analytic Poincare 


Key words and phrases. Spiral, box dimension of unbounded sets, Minkowski content, 
planar vector field, Hopf-Takens bifurcation at infinity. 
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map associated to th e spiral trajectories near sing ular points and periodic 


Zubrinic &: Zupanovic(2008) 


orbits can be seen in 

The multiplicity of the Poincare map is related to the notion of cyclicity, 
that is the number of limit cycles that can be born after a small perturbation 
of the system. Furthermore, the Poincare map of a planar vector field gener¬ 
ates a 1-dimensional discrete dynamical system. In [Elezovic et a/. (2007)] it 
was shown that the box dimension of the 1-dimensional trajectory is related 
to this discrete dynamical system. There the box dimension of the corre¬ 
sponding orbit is studied for the classical saddle-node and period doubling 
bifurcations. Using these results, classical theorems about these bifurcations 
were extended. Further extensions to 2-dimensional discrete dynamical sys¬ 
tems and applications to continuous dynamical systems were obtained in 
Dmitrovic Horvat(2012)| . 


It is known that limit cycles can also be generated from a polycycle, which 
is an ordered collection of singular points (vertices) and bi-asymptotic trajec¬ 
tories (edges) connecting them in a specified order. Remark that an isolated 
singular point is a special case of a polycycle. The next simplest case is 
a saddle-loop, that is a polycycle with only one vertex and one edge. The 
Poincare map near a saddle loop, although it is not analytic, shows its cyclyc- 
ity (see |Roussarie(I998)| , [Zhao Sz Wang(2009)| ). In [Mardesic et a/.(2011)] 
this was investigated from the point of view of fractal geometry. The classi¬ 
cal box dimension was not fine enough to distinguish between all the cases 
which could appear, so a generalisation called the critical Minkowski order 
has been introduced. 

As limit cycles can also be born from a point or a polycycle at infinity, it 
makes sense to generalise the previous results to this case. It is also inter¬ 
esting to study the problem of Hopf-Takens bifurcation of polynomial vector 
fields at infinity from the fractal point of view. Related problems have been 
studied in [Caubergh et a/. (2011)] , [Blows et a?. (1993)] and |Gine(2001)j . 

In this paper we deal with vector fields possessing spiral trajectories tend¬ 
ing to infinity. The provided visualisations clearly show that in the case of 
a weak focus at infinity such trajectories exhibit an almost “planar” na¬ 
ture. We measure this phenomenon using the box dimension of trajectories. 
Since the trajectories tending to infinity are unbounded, we have adapted 
the definition of box dimension to this case, since the usual box dimension 
is defined for bounded sets only. We do this using the geometric inversion, 
see Definition 12.31 below. 

Let us recall the definitions of Minkowski content and box dimension. By 
|D| we denote the n-dimensional Lebesgue measure of an open subset D of 
M"". Let A be a nonempty bounded subset of M", the e-neighbour hood of 
A in the Euclidean metric and s > 0. The upper s-dimensional Minkowski 
content of A is defined by 


M*^{A) = limsup —^ 
£^o e” * 


(1) 
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and we define analogously the lower s-dimensional Minkowski content of A, 
denoted by The upper box dimension of A is defined by 

( 2 ) 


AiuibA = inf{s > 0 : M.*^{A) = 0}; 


it is easy to see that we also have 


( 3 ) 


dims^d = sup{s > 0 : M*^{A) = oo}. 


The lower box dimension of A, denoted by dim ^^d. is defined analogously, 
using Ml{A) instead of M*^{A) in Q (and in Q). If both dimensions 
dim^^d and dim ^^d are equal, the common value is denoted by dim^A, 
and is called the box dimension of A (also known as Minkowski-Bouligand 
dimension, or limit capacity). If there exists d> 0 such that 0 < Ai'^{A) < 
JV[*^[A) < oo, we say that A is Minkowski nondegenerate, and Minkowski 
degenerate otherwise. (Note that if A is nondegenerate, it then follows from 
Q-Q and their counterpart for M-l{A) that dims^d exists and is equal 
to d.) If M'^{A) = M*'^{A), the common value is denoted by M'^{A), and 
called the Minkowski content. If moreover A4‘^{A) G (0, oo), then A is said 
to be Minkowski measura ble. For more information abo ut these notions and 
their generalisations see Zubrinic &: Zupanovic(2005)j , [Basic et a/.(2011)] , 
|Falconer(1990)| and [Mardesic et a/. (2011)] . 

In the sequel we will use the following notation. If /, g : M —>■ (0, oo) 
are two functions such that f{t) —)• 0 and g{t) —)• 0 as t —)• to {to can 
be oo as well), we write f{t) ~ g{t) as t —)• to if hmi_^tQ = 1. We 
write /(t) ~ g{t) and say that / and g are comparable as t —>• to if there 
exist positive constants ci ^2 such that ci 5 f(t) < /(t) < C 2 g{t) for all t in 
a neighbourhood of to- A function / : B —)• M”, V C is said to be 
Lipschitzian if \ f{a) — f{b)\ ~ |a — 6| for all a,b G V. 


2. Box DIMENSION AND MINKOWSKI CONTENT OF UNBOUNDED SETS 

2.1. Definition of box dimension of unbounded sets by geometric 
inversion. We start with a polynomial system 

(4) X = P{x) 

defined on M"'. Applying the change of variables u = x/|xp (here |x| is 
Euclidean norm, |xp = + • • • + x^), which is the well known geometric 

inversion of M” \ {0} with respect to the origin, after a short computation 
we arrive at the following system: 

(5) ii = \u\‘^P{u) — 2u{u ■ P{u)), 
defined on M” \ {0}, where 

The geometric inversion is clearly involutive, so that x = u/|up. The 
right-hand side of ([^ is not necessarily a polynomial field in dependence 
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of ui,... ,Un, where u = (ui ,... ,Un)- However, note that the largest expo¬ 
nent of \u\~‘^ appearing within the component functions of P{u) is equal to 
k = degP := max, deg Pi. Hence, 

(6) ii = (^\u\‘^P{u) — 2u{u ■ P{u))j , 

is a polynomial vector field. 

For any set A C MA \ {0} we can define its geometric inverse with respect 
to the origin by 

$(H) = {^>(x) : X G H}, 

where <h(x) = As we have said, the mapping is involutive: = id. 

1^1 

If we denote the phase portrait ofQby'P = {rj:zG/} (the family of 
trajectories Fj), it will be convenient to define <h('P) by 

(7) HP) = {cl>(F0 :iGl}. 

It is clear that <h('P) is the phase portrait of ([^ on M” \ {0}. Hence, we have 
proved the following result. 



Figure 1. The bounded spiral r = (left) and the un¬ 
bounded spiral r = (right) both have the same box 

dimension equal to 8/5. The nucleus of the bounded spiral 
is at the origin whereas the nucleus of the unbounded spiral 
is at infinity. 
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Lemma 2.1. Let a polynomial vector field P in M"' be given, and let V he 
the phase portrait of Then there exists an explicit polynomial vector 
field in M"’, given by such that its phase portrait is equal to ^{V). 

In particular, if {Ci : i £ 1} is the collection of all limit cycles of a 
polynomial vector field, then there exists a polynomial vector field in 
such that {<^>((7^) : i £ 1} is the collection of all limit cycles of the new 
vector field. 

As we see, if C is a limit cycle of a polynomial system, then its geometric 
inverse <h(C') is also a limit cycle of a polynomial system. 

Definition 2.2. We say that the infinite point is a weak focus of a dynamical 
system in M"' if the origin is a weak focus of the system obtained by its 
geometric inversion. 

It will be convenient to extend the definition of the box dimension from 
bounded sets in M” to the case of unbounded sets. 

Definition 2.3. Let A be an unbounded set in M”, which is away from the 
origin, that is, d(A, {0}) = inf{|a| : a G A} > 0. Then clearly $(A) is 
bounded, and we define the upper box dimension of A by 

dim^A = dim5<h(A). 

Analogously for the lower box dimension. If both the upper and lower box 
dimensions of A coincide, we call it just the box dimension of A, and denote 
it by dims A. 

The definition of the box dimension of A does not depend on the choice 
of the origin. 

Proposition 2.4. Assume that A is a given subset ofMP, and 0 ^ A. Assume 
that also w ^ A. Let $ be the geometric inversion with respect to the origin, 
and T the geometric reflection with respect to the point w £ M"', that is, 
^( 3 ^) = \x-w\2 • / = T o <I> : <h(A) —>■ 'I'(A) is a bi-Lipschitz mapping. 

In particular, see |Falconer(1990)| , we have dim 5 <I>(A) = dimB'L(A), and 
similarly for the lower box dimension. 

Proof. Let us first show that / is Lipschitzian. It suffices to show that 

(8) sup ||/'(x)|| < 00 . 

xS'I'jA) 

As the matrix norm || • || we take any operator norm, say oo-norm. First, 
by direct computation we see that 
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where I is the identity matrix, and x ® x = x ■ , with x understood as a 

column vector. Now, denoting a = <l>(x) we have: 

f'{x) = d/'($(x)) • d>'(x) 


a — w\‘^I — 2{a — w) ® {a — w) l®l 



Therefore ||/'(x)|| < C , where C is a positive constant and a = d>(x). 

Since a ^ A and —)■ 1 as a —)■ oo, the expression is bounded by 

a constant independent of a. This proves ®. 

On the other hand, f~^ = = ^oT, and we can show 

in the similar way that f~^ is Lipschitzian. Hence, / is bi-Lipschitzian. □ 


We can define the upper and lower s-dimensional Minkowski contents of 
A as the corresponding upper and lower Minkowski contents of <h(H), s > 0. 
We say that A is Minkowski nondegenerate (Minkowski measurable) if ‘h(H) 
is nondegenerate (Minkowski measurable). 


Remark 2.5. It is easy to get rid of the condition for A to be away from 
the origin. Indeed, if A is any set in MA, we can proceed as follows. Define 
Ai = An Bi{0) and A 2 = H \ Hi, and define 

dim^H = maxIdimsHi, dimB<I>(H2)}. 

It is easy to see that the upper box dimension so defined for unbounded sets 
satisfies the property of monotonicity (indeed, \i An B then ‘h(H) C ^{B), 
hence dim^ A = dim^ ‘^(H) < dim^ <h(H) = dim^ B), and the property of 
finite stability. See |Falconer(I990)| . 

Another basic property, as expected, is that the box dimension is pre¬ 
served for unbounded sets with positive distance from the origin that are 
bi-Lipschitz equivalent. 

Theorem 2.6. Let / : Vi —)• V 2 be a bi-Lipschitz map, where 0 ^ Hi and 
0 ^ ¥ 2 - If A C Hi, then dim^H = dimBf{A), and analogously, dim pH = 
dimp/(H). 


Remark 2.7. It is po ssible to constru ct a set H in M"' such that dim^H = n 
and dim ^H = 0, see |Zubrinic(2006) . 


Theorem 2.6 follows immediately from the following proposition, the proof 
of which we postpone. 


Proposition 2.8. Let Vi ^2 be two neighbourhoods of 00 in M"' such that 
0 ^ Hi and 0 ^ H 2 . The mapping / : Hi —?• H 2 is bi-Lipschitzian if and only 
if the mapping g : <h(Hi) —)• <I>(H 2 ) defined by 

(9) g{x) = o f o d>){x) 
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is bi-Lipschitzian. 

Proof of Theorem \2.6[ We have that C and using Proposi¬ 
tion 2.8 we obtain that ~ g'(<h(^)), with g defined by Q. Hence, 

dim^jd = dims<f>(74) = dims(7(<f>(74)) = dimB(<h o / o 

= dimB^>(/(^)) = dimB/(^). 

In the last equality we exploited the property of bi-Lipschitz invariance of 
the upper box dimension for bounded sets. Analogously for the lower box 
dimension. □ 


To prove Proposition 2.8, we start with the following two lemmas. 

Lemma 2.9. For any a,b € MF \ {0} we have 

|a — b\ 


ma)-m\ = 


|a| \b\ ’ 


where I I is the Euclidean norm. 


Lemma 2.10. Let f : Vi ^ V 2 be a bi-Lipschitz map, where 0 ^ Hi and 
0 ^ ¥ 2 - Then there exist two positive constants Ci and C 2 such that for all 
a € Vi, 

Cl < < C 2 . 

\a\ 

We omit the proofs of the previous lemmas and proceed to prove the 
proposition. 


Proof of Proposition 2^. Assume that / is bi-Lipschitzian, i.e. \f{a)—f{b)\ ~ 
|o — 6| for all a,b G Vi. Let x = <h(a) and y = <h(6) be any two elements 
from <l>(Vi). Using Lemma 2.9 we have 


|/(^>(x)) - /(^>(y))l |$(x) - ^{y)\ = 

Therefore, 

\fmx))-my))\ ^ 

\x — y\ 


x-y\ 


= \x-y\ |^>(x)| |^(y)|. 


kl \y\ 

|$(x)| my)\ 


\fmx))\\fmy))\ "'|/(ch(x))||/(ci>(y))|- 

Applying Lemma |2.9| on the left-hand side, and Lemma |2.10| on the right- 
hand side, we obtain 

|^>(/($(x))) - $(/(^>(y)))l ~ |x - y|, 

i.e. \g{x) — g{y)\ cz \x — y\ for all x,y G <h(Vi). 

The proof of the converse implication is similar, and therefore we omit 
it. □ 


If r is a smooth curve (typically, an unbounded spiral) in M” converging 
to infinity, which does not pass through the origin, then we can define its 
Minkowski content as follows. Assume that d = dim^r is well defined. 
Then we define Al'^(r) = Al‘^(<I>(r)). Note that the right-hand side is well 
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defined, since the set ‘h(r) is bounded. Furthermore, if we remove from 
r a portion of finite length, then the remaining part F^ has the same d- 
dimensional Minkowski content as F. This is due to the excision lemma, see 


Zubrinic(2006), Lemma 5.6]. 


Example 2.11. Let a be a given positive real number, and A = {k^ : k €N}. 
It is well known that the box dimension of d’(j4) = {k~°^ : k G N} is equal 
to 1/(1 + a), see e.g. [Lapidus &: van Frankenhuysen(2006)| . Therefore, 

dims A = —. 

1 “h Cfc 


The above example is a special case of the following result dealing with 
monotone strings C = {Ij) of infinite length, i.e. sequences of positive real 
numbers such that h ~ °° ih) nonincreasing. We do not 

require that /j —)• 0 as j —>■ oo. This string is associated with an unbounded 
sequence of real numbers A = (a^) defined by h- Conversely, it 

is clear that a nondecreasing, unbounded sequence of real numbers A = (a^) 
defines the string L = where Ij = a^+i — aj, and the string C = (Ij) 

is monotone if we require that Ij is nonincreasing. Note that here the set 
<I>(A) = {a'^^ : /c G N} is bounded, so that the classical box dimension makes 
sense. If we denote 


(10) flk — 0,p, Ofc+iJ ^ — ik'k) 

then 


(11) dims<l>(74) = dim^T' := inf/y > 0 : < co}- 


See Lapidus Sz van Frankenhuysen(2006)], where the right-hand side of (11) 


is taken as the definition of the upper box dimension of a general bounded, 
monotone string C = (^fc), denoted by D^' in this reference. Note that since 
< oo, then dims‘h(^) < 1. The following simple lemma provides a 
sufficient condition for a string associated with the geometric inverse of an 
unbounded set to be monotone. 


Lemma 2.12. Let A = (ak) be an unbounded, monotonically nondecreasing 
sequence of positive numbers. The string Cf = (fik), defined by (10), is 
monotone if and only if for each k > 1, 


( 12 ) 


Qfc+i Qfc+i ^ 2 

Ofc ak+2 ~ 


Furthermore, 

(AuibA = dim^Th 


Proof. It is easy to check that p.k+i < hk is equivalent with (12). 


□ 
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Example 2.13. For where a is a fixed positive number, the condition 

(12) is fulfilled, since 


Jk+l 


+ : 


= (l + fc-i)“ + (l + (fc + l)-')-“ > (l + k-T + C^ + k-T" > 2, 


\ —a 


— 1 \Q 




^k+2 


-1 


> 


where the last inequality follows from the elementary inequality t + t 
2 for t > 1. Therefore the conclusion of Example 2.11 is a special case 
of Lemma 2.12 since for a-strings C! = — [k + 1)~“)a:>i we have 

dim^T' = 1/(1 + a), see [Lapidus &: van Frankenhuysen(2006)| . 

Example 2.14. Let F be a spiral defined in polar coordinates by r = 
where > </?o > 0, and a is a given positive constant. Then <h(r) is an 
unbounded spiral defined by r = (/?“, where (f > (po > 0- We have 

2 

dims $(F) = max{l, ——}, 

1 + a 

see [Tricot(1995) p. 121]. Note that the nucleus of the spiral F is concen¬ 
trated near the origin, so that the nucleus of 4>(r) is concentrated at infinity. 
For a strict definition of the nucleus see [Tricot(1995)| , intuitively, it is the 
part where the e-neighbourhood of the spiral selfintersects. See Figure 

Example 2.15. Let a and [3 be two given positive constants. Let A be 
defined as the union of two spirals Fi and r 2 , defined in polar coordinates 
as follows: Fi ... r = (p~^ when p > 1 (bounded spiral tending to the origin), 
while r 2 ... r = when p > 1 (unbounded spiral, away from the origin). 
It is easy to see, using finite stability of the box dimension, that 

f 2 

dims A = max 1, —-——- 

[ 1 -|- minja, p| 

Starting from x = P{x), see (Q, using geometric inversion we arrived at 
u = P*{u) where 

(13) P*{u) = \u\‘^P{u)-2u{u-P{u)). 

It is clear that P** = P for each vector field P, since the geometric inversion 
with respect to the origin is involutive. It is easy to see that 

* : C^(M" \ {0},M”) ^ \ {0},M^) 

is a linear operator with real coefficients: for any A, p G M and F,G € 
Ci(M’^ \ {0},M’^) we have (AT + pG)* = XF* + pG*. 

Remark 2.16. If in (Q P{x) is a rational function (that is, the component 
functions are rational functions of Xj, j = 1,... ,n), then from (13) we see 
that P* (u) is also a rational function. The phase portrait of the system is 
the same (outside the origin) as for the polynomial system corresponding to 
d(x)P(x), where d(x) is the common denominator of all Pj(x}. Analogously 
for the system (j^. 

The following lemma deals with a special class of right-hand sides P(x) 
of a for which P*(u) can be easily computed. 
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Lemma 2.17. Let us consider the system ^ with P{x) = Rx — 'yxg{\x\), 
X G where 'j is a real constant and g : (0, oo) —)• M a continuous function, 
and R is an n X n real antisymmetric matrix: R~^ = —R. Then 

P*{u) = Ru + 'yug{\u\~^). 


Proof. The matrix R is antisymmetric if and only if Rx ■ x 
The claim follows from (13) and Pu = 
computation. 


u\ ^{Ru —'yug{\u\ ^)) 


= 0 for all X. 
after a short 
□ 


Example 2.18. In particular, if P{x) = Rx, where R is a, real antisymmetric 
matrix, then P*(u) = Ru, that is, P = P*. If P{x) = cx, where c is a real 
constant, then P*{u) = —cu. 


A typical example of a real matrix R satisfying the condition Rx ■ x = 0 
for all X G M"' is any diagonal block matrix containing either matrices of the 
form 


on the diagonal (here Xj G M), or zeros. 

If we deal with an ODE in the complex phase space C"": 

z = P{z), 

where z = {zi,..., ZnY, then introducing the new variable u = zl\z^\, where 
Izp = obtain 

u = P*{u), 

where 

P*{u) = \u\^P{u) - 2u Re(n | P{u)). 

Here we define {u\v) = X]j=i 

Lemma 2.19. Let P{z) = Rz — 'yzg{\z\), where R = diag(iAi,..., 

Xj G M, 7 zs a given complex number and g : (0, oo) oo is a continuous 
function. Then 

P*{u) = Ru + 7 zzg'(|zz|“^). 

Proof. Note that {Ru\u) = iYlj ^ that Re(i2zt|zz) = 0. The rest 

of the proof is the same as in the proof of Lemma |2.17[ □ 


We shall also often need the following technical lemma, dealing with pla¬ 
nar systems of weak focus type near the origin; 


(14) 


X = -y+p{x,y) 
y= x + q{x,y). 


A typical situation is when p and q are analytic functions with McLaurin 
series containing quadratic or higher order terms only. It is an extension of 
Lemma 12.171 in the case of n = 2. 
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Lemma 2.20. 

equal to 


The system obtained from (If) by geometric inversion is 


(15) 


ii = —V + {v‘^ — u^)p — 2uv q 
V = u + {u^ — v‘^)q — 2uvp, 


where p = p{^^, and q = q{^^, 

The corresponding general result concerning the box dimension of spiral 


trajectories of (15) can be seen in Theorem 3.3 


2.2. Definition of box dimension of nnbounded sets in the plane us¬ 
ing the Riemann sphere. Assume that A C is an unbounded set such 
that the origin is not its accumulation point, so that d>(A) is bounded. Let 
S be the Riemann sphere in of radius 1/2 with the center at (0,0,1/2). 
Let TTs : ^ S he the stereographic projection of the plane to the sphere. 

The following result shows that the box dimension of A (defined via geo¬ 
metric inversion) is the same as the the box dimension of its stereographic 
projection. Here the set 7r5(A) contained in the Riemann sphere, is viewed 
as a subset of M^, i.e. its box dimension is computed via its e-neighbourhood 
in 


Proposition 2.21. Let A be an unbounded set in such that 0 ^ A. Then 
dim^A = dimB7r5(A), analogously for the lower box dimension. Further¬ 
more, if A is Minkowski nondegenerate, so is 7r5(A). 

Remark 2.22. If A' is a subset of the Riemann sphere S, then it is possible 
to dehne its e-neighbourhood on the manifold S with respect to its metric 
and the resulting surface measure. This permits us to define the new (upper 
or lower) box dimension in the usual way. It can be shown that the box 
dimension of A' with respect to the manifold S is the same as the box 
dimension of A' with respect to i.e. with respect to its e-neighbourhood 
in 


The proof of Proposition 2.21 rests on the following lemma. 


Lemma 2.23. Let Hro(O) the disk of radius tq G (0,1) in the plane, 
and F : (0) —)• S defined as composition of inversion and stereographic 

projection, i.e. F{x) = for x h, and F{0) = N (the north pole of 

S). Then the mapping F is bi-Lipschitzian. 


Proof. Let P : S' —)> be the orthogonal projection. It is easy to see 

that V = P(Rro(0)) is a neighbourhood of the north pole N = (0,0,1) 
on the Riemann sphere, and strictly above the equator, due to tq < 1. 
It is then clear that the restriction of the projection P\v is a bi-Lipschitz 
mapping. Therefore, it suffices to show that G = P o F is bi-Lipschitzian, 
since then F = P~^ o G will be bi-Lipschitzian as a composition of bi- 
Lipschitz functions. 
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Figure 2. The unbounded spiral f{(p) = projected to 
the Riemann sphere of radius 1/2. 


We pass to polar coordinates (r, cp) in the plane. It is clear that G has 
the form G{r,(p) = {g{r),(p). Fixing any vertical semi-plane defined by 
given if, it intersects the Riemann sphere in a semi-circle above the r-axis. 
For a given r > 0, the strait line joining r~^ with N interesects the circle 
in the point the horizontal component of which is equal to g{r). An easy 
calculation shows that 


(16) 


g{r) 


r 

_|_ 1 


Since g'{r) = we have that ^'(O) = 1 and g'{ro) > 0, hence g'{r) G 

{g'{ro), 1). This shows that G : V ^ G{V) is bi-Lipschitzian. □ 


Remark 2.24. It is easy to see that the bound rg < 1 in Lemma 2.23 
optimal, since g'{l) = 0. 


IS 


" be a given spiral A in the plane, ip > (po > 0. 

is a spiral defined by 


Example 2.25. Let r = p 
Then G{A), with G from the proof of Lemma 2.23 

P = 


— y 




In particular. 
The following 


and both of them are bi-Lipschitz equivalent. 

they both have the box dimension equal to d = maxjl, 
proposition shows that they both have t he same d-dimensional Minko wski 
content (its value has been computed in Zubrinic & Zupanovic(2005) , see 
below). 


Proposition 2.26. Assume that a G (0,1) is a given constant. Let r = f{p) 
he a spiral such that f{p) Cp-^, f{p) —aGp " ^ as p ^ oo, and 
there exists M > 0 such that \f”{p)\ < Mp~°‘, for all p > po > 0. Define 
the new spiral p = ■ Then both spirals have the same box dimension 

d = 2/(1 -|- a), and the same d-dimensional Minkowski contents. 
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The proof of this propositio n is a direct consequence of 


Zubrinic & Zupanovic(2005), Theorem 6] that we state here in a simplified, 


but equivalent form (stated also in Korkut et a/. (2009), Theorem 3]). 

Theorem 2.27. (Minkowski measurable spirals) Assume that f : [(/Ji, oo) —)> 
(0, oo) is a decreasing -function, and tpi > 0. Assume that there exists 
the limit 


m = lim 




— Ol\1 ’ 


and m > 0. Let there he a positive constant M such that \f'\ip)\ < 
for all if > ^ 1 - Let T be the graph of the spiral p = f{ip) with a G (0,1) and 
d = 2/(1 + a). Then dim^T = d, the spiral is Minkowski measurable, and 
moreover, 

1 + 0 


(17) 


M'i(r) = m'^7r(7ra)-2“/(i+“) 


1 — 0 


Now the proof of proposition |2.26| follows by showing t 
^(+) = the same properties as f{(f) in Theorem 

with a different value of M >0. 


rat t he functon 
2.271 eventually 


Remark 2.28. The analogous construction as in Proposition 2.21 and in 
can be performed starting with ML instead of and using 


2.23 


Lemma 

the Riemann sphere S'"' in of radius 1/2, centered at (0,..., 0,1/2). It 
suffices to use the stereographic projection vr^n : M" —)• S". 

Example 2.29. Let us suppose that the Riemann sphere is of radius R. 
Proposition 2.21 still holds which is easy to see using analogous arguments 
as in Lemma 2.23 In this case we take rg G (0,1/(2R)) and get 

, , m^r 
~ 4RV + 1’ 

On the other hand, the Minkowski content of the new spiral p will be afected 
with the radius of the Riemann sphere. Concretely: 

Let be the Poincare sphere in of radius R, i.e. 5^ = {{X, Y, Z) G 

We shall project onto the sphere from the 
(x, y) plane placed tangentially at the north pole. We are interested how 
the box dimension of a focus type spiral is affected with geometric inversion 
and projection onto the Poincare sphere. 

Proposition 2.30. Let Pi ... r = /(y?) be a spiral of focus type such that 

f{g>)c.g>--, |/'(^)|c.^-“-\ |/"(^)|<M^-“ 

as ip ^ oo, for some positive constants a and M. Firstly we geometrically 
invert this spiral, and then project it on S'^. We shall denote this new spiral 
in with + 2 . Then 

2 + 0 


(18) 


dims+2 = 


1 + 0 
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Proof. In cylindrical coordinates in the above described map of the spiral 
Fi is given by 




R 








We will use 


Zubrinic &: Zupanovic(2005)[ Theorem 5(b)]. In our case we 
have r = R — F(ip), where 


F{ip) = R- 


R 


ViTWfW 

and Z = y/R‘^ — = g{\R — r\), for g{t) = V 2Rt — if. 


Now we have to check the conditions from Zubrinic & Zupanovic(2005) 
Theorem 5(b)]: 


g{t) ~ g'{t) 




We can see that g meets the conditions for (3 = 1/2. Now we check the 
conditions on F\ 

F{p) ~ 

\F\^)\ = 


(1 + i ? 2 /(^) 2)^/2 


R^ 


-2a 


^ (l + ^2j2)5/2 d/'l' + f\f"\ + + 2R^f\f?) < 

for ip sufficiently large. So w e can see that F meets the conditions of the 


Zubrinic & Zupanovic(2005), Theorem 5(b)] for a 2 = 2q: and we have the 

2 + 02/3 2 + a 


conclusion 


dims r2 = 


l + a2(3 l + a' 


□ 


Remark 2.31. For a > 1 the box dimension of the spiral Fi is one. This 
shows that projecting a spiral of trivial box dimension to the Poincare sphere 
we can get a spiral of nontrivial box dimension. 


Remark 2.32. As the Poincare sphere is usually represented by orthogonally 
projecting it on the xy-plane, this will further affect the box dimension of a 
spiral defined in Proposition 2.30[ After projecting it, we will have a limit 
cycle type spiral F 3 ...r = R — F (ip) in the plane and its box dimension will 
be reduced to 


2 + 2a 
l + 2a' 


dims F 3 
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Figure 3. The unbounded spiral f{(p) = projected to 
the Poincare sphere of radius 1. 


For a focus type spiral Fi of nontrivial box dimension i.e. for a G (0,1) 
we have the next relations between box dimensions: 

dims r2 = 1 + ^ dims Fi 


dims Fs 


3 —dimsr2’ 4 —dimsFi" 

3. Weak focus at infinity 


3.1. Weakly damped oscillator. Let us consider a weakly damped oscil¬ 
lator 

y+ C'y"(y)^ + y = o, 

where a is even and /3 odd positive integer, and C is any positive constant. 
It is well known that it is globally stable. It is equivalent to the following 
planar system: 

X = —y — Cx^y°‘ 


y = 


X. 


All nontrivial trajectories F, corresponding to t > 0, are spirals converging 


clockwise to the origin, and the origin is the weak focus. Using Lemma 2.20 
we conclude that the corresponding system obtained by geometric inversion 


IS 


ii = —V + C 


(19) 


-I- 


V = u + 2C 


^0+lya+l 

(n2 +u2)«+/3- 


All trajectories corresponding to t > 0, starting outside the origin, are of the 
form ‘^(F) for some F as above. The spirals <I’(r) are converging clockwise 
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to infinity, and the infinity is the weak focus. The system (19) becomes 
polynomial after multiplying the right-hand sides by (u^ + 

u = -v(u^ + + Cu^u^fu^ - v^) 

1201 


Theorem 3.1. Any nontrivial trajectory T of the system (20), corresponding 
to t > 0, is a spiral converging to infinity, and 

dim^T = 2(1 - 

y q; + 

Furthermore, the spirals are Minkowski nondegenerate. 


This follows immediately from Pasic et a/. (2009), Theorem 7]. 

3.2. Lienard systems. Let us consider the following Lienard system: 


N 


N 


x = -y + '^ a2iX^'’ + ^ a2j+ix^*+^ 
i=l i=k 

y= X. 

Here we assume that a 2 A:-i-i 7 ^ 0, which means that this is the first nontrivial 
coefficient on the right-hand side having odd index. The system obtained 
by geometric inversion is 

ii = —V + — u^)p(u) 

V = u — 2uvp{u), 

where 


( 21 ) 


N 


V{u) = '^0,2i- 


u 


2i 


N 


if + X]“2i-ri- 


,2i+l 


i=l ^ i=k 

Multiplying the right-hand sides of ( |2l[ ) by {vf + 0 "^)^, the system becomes 
polynomial, retaining t he same phase portrait outsid e the origin. An im¬ 
mediate consequence of Zubrinic &: Zupanovic(2008), Theorem 6 ] is the fol¬ 
lowing result. 


Theorem 3.2. If in (21) we have a 2 k+i / 0, then for any initial point 
{uo,vo) 7 ^ ( 0 , 0 ) we have that the corresponding trajectory T = r(uo,uo) 
starting from that point is a spiral tending to infinity, and 

dims r = 2(1 — 


2k+ 1 

Furthermore, T is Minkowski nondegenerate. 

As we see, unbounded spiral trajectories (in fact, semitrajectories, i.e. 
starting from initial point) of Lienard systems can achieve box dimensions 
with values from the following set only: 

r A 8 12 16 20 . 

^ ^2k + l ^ ^3’ 5’ 7 ’ 9 ’ 11’ 
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For analytic systems, these are the only values of box dimensions that un¬ 


bounded spiral trajctories can achieve, see Zubrinic & Zupanovic(2008) . 

A more general result can be stated in terms of the Poincare map at 
infinity. We deal with the system (14) such that p{x,y) and q{x,y) are 
analytic functions of the form 


( 22 ) 


Pix, y) = "^Phix, y), q{x, y) = '^ qk{x, y) 


k=2 


k=2 


where pk and qk are homogeneous polynomials of k-th. degree. The Lyapunov 
coefficient of a system at infinity is defined as the Lyapunov coefficient at 
the origin of the system obtained by geometric inversion. The Lyapunov 
coefficient near the weak focus is defined as the coefficient of the leading 
term of the Taylor expansion of the displacement function. T he following 
result follows immediately from Zubrinic &: Zupanovic(2008), Theorem 6]. 


Theorem 3.3. Let T be an unbounded spiral trajectory, away from the ori¬ 


gin, associated to the system (15). Assume that p{x,y) and q{x,y) are ana¬ 


lytic functions as in (22). IfV 2 k+i is the first nonzero Lyapunov coefficient 


of (15) at infinity, then 


dims T = 2(1 — 


1 


2k + I 

Furthermore, T is Minkowski nondegenerate. 


3.3. Classical Hopf bifurcation. The classical Hopf bifurcation is dehned 
by the following system for k = 1: 


(23) 


X = -y - x (^{x‘^ + y^)^ + 

y= 2: - y (^(x^ + y^)^ + , 


where a is the bifurcation parameter. The corresponding spirals T are con¬ 
verging clockwise to the origin. Using Lemma 2.17 (here R is the symplectic 
2x2 matrix and y(r) = r"^^ + o) we have that the related system obtained 
from (23) by geometric inversion is 


(24) 


ii = —V + u ^ 

V = u-^v -j- v^)~^ -b aj , 


and the corresponding spirals are converging clockwise to infinity, which is 
the weak focus. The corresponding polynomial system 


(25) 


u = —v(u^ -b -b u ^1 -b a(u^ -b 

y= uiu"^ + + V (1 + -\-, 


has the same phase portrait as (24) outside the origin. 
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(26) 


In polar coordinates (r, (/?) system ( |23[ ) has the form 

r = —r(r^^ + a) 
ip=l. 


For a < 0 the limit cycle is born off the origin, r = (—while system 
(24) in polar coordinates {p,^p) has the form 

(27) P = + 

ip=l. 

In this case, for a < 0 the limit cycle is born off infinity. Here r = (— 
and r —?• oo as a —7- 0_. The system (27) is clearly the same as the one ob- 

cos 1/3 


_ Sinip rp 

y p ■ 


tained from (23) by introducing the coordinates (p, ip) dehned via x = 


re following result shows that the box dimension ‘recognizes’ 
the Hopf bifurcation. 


Theorem 3.4. Let a = 0 in the bifurcation problem (24) or (25). Then any 
unbounded spiral trajectory T, away of the origin, has the box dimension 
equal to 

dimBr=^, 

2/t + l’ 

and is Minkowski measurable. For all the other values of a the box dimension 
is trivial, i.e. equal to 1. 


This is an immediate consequence of Zubrinic &: Zupanovic(2005), The¬ 
orem 7]. 


4. Hopf-Takens bifurcation at infinity 


Using geometric inversion and results from [Zubrinic &: Zupanovic(2005) , 
chapter 4, we shall study fractal properties of Hopf-Takens bifurcation oc¬ 
curring at inhnity. For a standard generic Hopf-Takens bifurcation we have 
the normal form; 



where (oq, •. •, Ui-i) £ is fixed. In sequel we will consider only, since 

the case is treated similarly. In case the normal form in polar 
coordinates is given by 


(28) 


r = r r 


21 


l-i 

+ E 

i=o 


a^r 


2i 


(p = l. 
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Geometric inversion with p := 1 /r yields a new system of differential equa¬ 
tions 


(29) 


P = 



ip = l. 

Now it is easy to see that t he following analogous versions of 


Zubrinic & Zupanovic(2005), Theorems 9 and 10 ] are valid. 


Theorem 4.1. (The case of focus). Let T he a part of a trajectory of (29) 
near infinity. 

(a) Assume that oq 0. Then the spiral T is of exponential type, that is, 

comparable with p = and hence dim^ T = 1 . 

(b) Let k be fixed, 1 < k < I, ao = ■ ■ ■ = Ok-i = 0, Ok 0. Then T is 

comparable with the spiral p = , and 

dimer =^. 

2A; + 1 

Theorem 4.2. (The case of limit cycle). Let the system (29) have a limit 
cycle p = a of multiplicity m, 1 < m < 1. By Ti and T 2 we denote the 
parts of two trajectories of (29) near the limit cycle from outside and inside 
respectively. Then the trajectories Ti and r 2 are comparable 

(a) with exponential spirals p = a± e~^‘^ of limit cycle type when m = 1, 
for some constants /3 7 ^ 0 (depending only on the coefficients Oi, 0 < i < 
l-l), 

(b) with power spirals p = ai: when m > 1. 

In both cases we have 

dims Tj = 2-, i = 1,2. 

m 


Remark 4.3. In Theorem 4.2 the parts of trajectories we are observing are 
contained in an open ring around the limit cycle which is a bounded set 
that does not contain the origin. As geometric inversion <1> is bi-Lipschitzian 
on such sets, Theorem |4.2| is a direct consequence of Theorem 10 from 


Zubrinic & Zupanovic(2005) 


Remark 4.4. From Theorem 4.2 we know that for (29) each spiral trajectory 
of limit cycle type has box dimension from the set 


T>i = {2 - — : m G N} = 
m 


3 5 7 9 
2’ 3’ 4’ 5’ 


See Zubrinic &: Zupanovic(2008), p. 958]. 

Let us have a look at the inversion of a standard Hopf bifurcation in polar 
coordinates, i.e. system (29) for Z = 1: 

P = + ao) 

ip = l. 


(30) 
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Viewing ao as a bifurcation parameter, we have the three following pos¬ 
sibilities. 

(1) For oo < 0 the trajectories of (30) are given with 


p{v) 


- —\/-ao 

an V 


C>0, (/?GMorC'<0, (/?< 


In(aoC) 

2ao 


We can see that we have a strong focus at infinity and the circle p = 
(—is the limit cycle for trajectories from inside and outside near 
the circle. The corresponding spirals near infinity are comparable with 
p = while the spirals near the circle are comparable with p = 

(—All these spiral trajectories are of exponential type and 
hence of box dimension equal to 1. See Figure]^ left. 

(2) For ao = 0 the trajectories of ( [30l ) are given with 

p{}p) = \/—2ip + C, C G M, (/?<—. 


and infinity is a weak focus with dim^ F = 4/3 where by F we denote a part 
of the trajectory near inhnity. See Figure]^ middle. 

(3) For ao > 0 the trajectories of ( fSO] ) are given with 

PiV’) = -J-ao + agC'e-2“o¥>, C > 0, y? < , 

ao » 2ao 


Inhnity is a strong focus and all the trajectories near inhnity are comparable 
with the spiral p = of exponential type, and hence have box dimension 

equal to 1. See Figure]^ right. 

Let us now consider the case I = 2 in (29): 


(31) 


p = -p{p -h ao aip 
ip = l. 



Figure 4. Trajectories of the system (30). Left; ao = — 1/25 
with a limit cycle born from inhnity; Middle: ao = 0 with 
weak focus at inhnity; Right: ao = 1/25 with exponential 
spiral at inhnity. 










FRACTAL ANALYSIS OF HOPF BIFURCATION AT INFINITY 


21 



Figure 5. Trajectories of the system (30) drawn on the 
Poincare disc. For oq = 0 (middle) the box dimension around 
the equator would be 1+^| = 5/3 but after projecting the 

half sphere onto the disc the box dimension is reduced to 
2 _ 3 

3 - 5/3 2 - 


Let US fix the value ai = —2 and consider ao as a bifurcation parameter. 
Since it is clearer to see what is happening, the phase portrets will be drawn 
on the Poincare disc. 





Figure 6. 



Figure 7. 
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(a) When ao < 0 all box dimensions are equal to 1 because all the trajec¬ 
tories are of exponential type, see Figure left. 

(b) For ao = 0 we have a weak focus at infinity and any part of a trajectory 
F near infinity has box dimension equal to d = 4/3 (power case), whereas the 
part near the limit cycle r = 1/\/2 has box dimension equal to 1 (exponential 
case), see Figure middle. Actually, because the process of projecting 
onto the Poincare disc affects the box dimension (see Remark [2(3^ , the 

trajectories on the figure near the equator have box dimension equal to 

4 _ 3 

4-d 2- 

(c) For ao G (0,1) we have two limit cycles of multiplicity one, and all 
box dimensions are eqnal to 1 (exponential case), see Figure]^ right. 

(d) For ao = 1 we have a limit cycle r = 1 of mnltiplicity two, and all 
trajectories near the limit cycle (either inside or outside) have box dimension 
eqnal to 3/2 (power case), see Figure]^ left. On the other hand, trajectories 
near the equator have box dimension eqnal to one (exponential case). 

(e) For ao > 1 box dimensions of all trajectories are equal to one (expo¬ 
nential case), see Figure]^ right. 
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